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Abstract. We consider the steady-state Navier-Stokes equation in the whole 
space driven by a forcing function /. The class of source functions / under 
consideration yield the existence of at least one solution with finite Dirichlet 
integral (||VL'^||2 < oo). Under the additional assumptions that / is absent 
of low modes and the ratio of / to viscosity is sufficiently small in a natural 
norm we construct solutions which have finite energy (finite norm). These 
solutions are unique among all solutions with finite energy and finite Dirichlet 
integral. The constructed solutions are also shown to be stable in the following 
sense: If U is such a solution then any viscous, incompressible fiow in the whole 
space, driven by / and starting with finite energy, will return to U . 



1. Introduction 

The classical theory of viscous, incompressible fluid flow is governed by the fa- 
mous Navier-Stokes equations: 

(1.1) ut+u-Vu + Vp^v/\u + f 

m(0) = Mo V • u = 

A large area of modern research is devoted to deducing qualitative properties of 
solutions for these equations when they are complemented with initial and boundary 
conditions and certain restraints are placed on / and v. The investigations in 
this subject are too numerous to attempt to list here so we will limit ourselves to 
discussion directly related to the topic of this paper: the steady state Navier-Stokes 
equation in the whole space 9? . 

A steady state (sometimes called stationary in the literature) solution U of the 
Navier-Stokes equation is one for which dtU — 0, that is the solution is constant 
in time. Such solutions solve the following PDE, which will be our main point of 
investigation. 

(1.2) U -VU + Vp = iyAU + f 

V • [/ = 

For our purposes this PDE is supplemented with the idea that U tends to zero as |x| 
becomes large, made precise by working in functions spaces which arc completions 
of smooth functions with compact support. 
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Roughly speaking, the investigation of steady state solutions can be broken into 
two regimes: bounded and unbounded domains. In the former situation much 
progress has been made using a Poincare type inequality (||C/|j2 < C||VC/||2) to 
deduce quickly that solutions have finite energy. In the case where there is no 
Poincare inequality it is desirable to to find conditions on / which will guarantee, 
a priori, finite energy of a solution. One of the benefits of establishing existence of 
solutions with finite energy in unbounded domains is that many of the techniques 
developed using the Poincare inequality can be applied, but the consideration of 
unbounded domains is not a needless complication. Indeed, many physical problems 
are best stated in the whole space or in exterior domains where there is no Poincare 
inequality. Moreover, the situation in the whole space is theoretically important as 
Leray observed in his seminal paper outlining modern analysis of the Navier-Stokes 
equations [2T: 

L 'absence de parois introduit certes quelques complications concer- 
nant I'allure a I'infini des fonctions inconnues, mats simplifie beau- 
coup I'expose et met mieux en lumigre les difficultes essentielles; 
The main goal of this paper is to develop a new technique which will allow, 
with certain conditions on /, the construction of solutions for the steady state 
Navier-Stokes equation in the whole space with finite energy. The assumptions we 
impose on / limit the amount of low frequency information and require that the 
ratio of / to I' is small in a natural norm. Once we have established the finite 
energy of solutions we deduce uniqueness in the class of solutions with finite energy 
and prove these solutions are stable in a strong sense referred to in the literature 
as nonlinearly stable. In other situations (not the whole space) non-uniqueness 
for solutions of (|1.2p has been demonstrated for solutions with / large compared 
with V (see [33], Chapter 2, and references therein) so we suspect the smallness 
assumption we make on / is necessary and natural for this result. It is currently 
unknown if the assumption on low frequency information is natural or a byproduct 
of our technique. 

1.1. Statement of Results. Modern analysis of the steady state Navier-Stokes 
equation in unbounded domains can be traced back to |20j and [26j , these ideas were 
further developed in [5j, [S], and |10| which work with the notion of a physically 
reasonable solution. The authors were concerned with the physically interesting 
problem of solutions in exterior domains of which the whole space is a special case. 
We summarize (perhaps too succinctly) these works with the following idea: If / 
is such that the Dirichlet integral ||VC/||2 is finite {\x\f e £^ is sufficient) then 
there exists a unique physically reasonable solution U in an exterior domain. This 
solution is physically reasonable in the sense that it approaches a constant (possibly 
non-zero) vector field like as becomes large and the uniqueness is among 

all such functions. These ideas were expanded further in [2], [3], [12], [13], [Hj, and 
[22] . The methods rely heavily upon analysis of the Green's function for the domain 
in question and are quite different from the approach presented in this paper. 

Our construction of solutions with finite energy is based on a well known formal 
observation: if $ is the fundamental solution for the heat equation then <^{t, ■) dt 
is the fundamental solution for Poisson's equation. Using this idea it is possible 
to make a time dependent PDE similar to the Navier-Stokes equation with / as 
initial data with a solution that can be formally integrated in time to find a solution 
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of ()1.2|) . At this point our analysis turns to the theory of energy decay for fluid 
equations. If the decay of the new time dependent PDE is fast enough (an integral 
over all time converges) we can deduce a finite energy bound for (|1.2p . It is known 
that the decay rate of solutions for parabolic PDEs in the whole space is intimately 
related to the shape of the initial data near the origin in Fourier space; the low 
frequency assumption we make on / is enough to guarantee the convergence of the 
required time integrals. This idea is further outlined in Section 2 and made precise 
in Section 3. 

A particularly useful technique for estimating energy decay is the Fourier Split- 
ting Method which was used in [31j to establish energy decay for initial data 
uq ^ and later for initial data uq G in [25]. Other works in this area 
include [2], ^7\, [ig, [22], [23], [24], [27|, [2g, [29], [30], [35], and ^36]. In essence 
we are trading bounds on the Green's function for energy decay theorems which we 
base on the Fourier Splitting Method. The assumption on low frequency informa- 
tion is stronger then the classical assumptions but our conclusion is stronger. Of 
course the previous results consider the more complicated cases of external domains 
which are not handled within but we hope that with decay theorems for external 
domains one can use the technique presented here to obtain similar results. We 
now state precisely the main theorem proved. In the following statement is 
the completion of smooth divergence free functions of compact support under the 
norm ||V • II2 and = n . Also, X = (H^)' n L^. The requirement f £ X 
implies the classical assumptions f £ L"^ and ||VC/||2 < 00. The later is known as 
a finite Dirichlet integral and is sometimes implied by the restriction \x\f G in 
the literature. 

Theorem 1.1. Let M > and f £ X satisfy the following assumption: 

(A) There exists a po such that /(^) — for almost every \^\ < po 
Then there exists a constant C{pq, M) so that if\\f\\x ^ C{po, v, M) the following 
hold: 

(i) The PDE has a weak solution U G H^. It is a weak solution in the 
sense that for any divergence free function of compact support (f), 

(1.3) <[/• VC/, (/)> < VC/, V(/. >=</,(/) > 

(ii) This solution satisfies \\U\\2 < M and ||VC/||2 < '^^"""II/Hat- 

(iii) This solution is unique among all solutions which have a finite norm 
and satisfy |jVC/||2 < ^^^"'"ll/llx- 

Remark 1.2. The behavior of the constant C{po, i^, M) allows large / when the 
Reynolds number is small (see remark r3.16|) . In this work we assume the Fourier 
transform of / is zero in some neighborhood of the origin, this corresponds to 
exponential decay for the heat flow starting with initial data /. It is possible to 
relax the hypothesis so the heat flow is algebraic but fast and not signiflcantly 
change any of the proofs presented here. This could be accomplished by requiring 
/ to behave like the polynomial |^|^ near the origin in Fourier space where p is some 
sufficiently large number (see pQ). 

When a solution is known to have finite energy the well developed energy stability 
arguments can be applied when / is small measured against see [6], [7], [32], and 
[33j . This is argued by showing solutions of a nonlinear parabolic PDE, similar to 
the Navier-Stokes equation and found by subtracting the steady state, tends zero as 
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time becomes large. Again we tm'n to energy decay methods, specifically a method 
applied to the Navier-Stokes equation in [25], to show the decay. The method 
consists of estimating the high and low frequencies of the solution separately and 
the estimates of the high frequency rely on the Fourier Splitting Method. Section 
4 is dedicated to establishing the following theorem. 

Theorem 1.3. Let f satisfy the assumptions of Theorem \L1\ and be such that 
||/||x is less then the constant given by the theorem. There exists another constant 
C(i^) such that \\f\\x ^ implies U is stable in the following sense: if Wq G L^ 

is a perturbation and u is a solution of the Navier-Stokes equation il.l]) with initial 
data Wq + U which satisfies, for any T > 0, 

ueL°° (0, T; i^) n £2(0^ T; Hi) 

then 

(i) For every e > there is a S > such that 

||wo||2 < S implies sup \\u{t) — U\\2 < e 
teR+ 

(ii) u{t) tends to U as time becomes large, that is 

lim \\u{t) ~ U\\2 = 



1.2. Notation. Unless otherwise noted all integrals in this paper are taken over the 
whole space E.^, C§° denotes the space of smooth functions with compact support. 

V = {0eCo°°|V-0 = O} 

L^ — { completion of V under the norm || • ||p} 

HI = { completion of V under the norm || • ||^i} 
iHly = {dual of Hi} 
II • \\x = max{|| • II2, II • ll(iji),} 

We will typically write an element / e {Hi)' as "/" when we really mean the map 
"(/) — >< f,(j) >" . To denote general constants we use C which may change from line 
to line. In certain cases we will write C{a) to emphasize the constants dependence 
on a. In a similar way we write p to denote general potentials (used to describe the 
pressure, one instance of p may not be the same as another even on neighboring 
lines). The variable ^ is reserved for working in Fourier Space. 



</,.?>= J f-9 

l|-||ffi=l|V-||2 

Hi = Lin Hi 
x = Lln{Hiy 

m = J f{x)e-^^"-<dx 
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2. Preliminaries 

Existence of weak solutions for ()1.2p is well known, see for example [5], [8], [TT] . 
[in], [in] J [20] J [31] • A typical approach to constructing weak solutions for this 
PDE is to construct approximations with the Galerkin Method and use a priori 
bounds with the Banach-Alaoglu theorem to find a subsequence of approximations 
converging weakly to a possible solution. Some stronger compactness property is 
then used to pass the sequence through the nonlinear term and establish the limit 
is indeed a solution. A good a priori bound for this approach, and a bound we will 
rely on throughout is: 

(2.1) m\\i<'^-'\\f\\x 

This is essentially the classical assumption that U has a finite Dirichlet integral but 
we derive it from our assumption f € X using the estimate 

|</,C/>|<||/||x||VC/||2 

The bound (|2.ip is proved formally by multiplying (|1.2p by U then integrating by 
parts. Noting the specific form of the nonlinearity, 

(2.2) <wv;7,;7>=o 

This relation holds when V • w = and the integral is absolutely summable, it 
can be proved for functions of compact support using integration by parts then 
extended to other classes of functions with a density argument. It holds in three 
dimensions when w G and U E H^, or when U E L^. O and w E since 
either assumption implies summability. 

Fix / and C/ as a solution to (|1.2p {U does not depend on time), we would like 
to find conditions on / which guarantee ||J7||2 < oo. One of the key steps of our 
approach is to establish "fast" decay of solutions to the system 

(2.3) Vt + U ■\'v + Vp = vAv 

v{0) = f V • w = 

Formally, if v is a solution of ()2.3p then U — v{t) dt solves 

U -VU + Vp = lyAU + f 
V • {/ = 

Recall we have fixed U earlier and it is also a solution for this PDE since it satisfies 
(|I.2p . As this PDE is linear and V • C/ = 0, solutions are unique and we may 
conclude U = U. Using Minkowski's Inequality for integrals we can see how the 
decay of v relates to the norm of U: 

/>00 /"OO 

||C/||2- II / V{t)dt\\2 < / \\vit)\\2dt 
Jo JO 

In summary, if ||i;(i)||2 < C(f + 1)^^^ with /? > 1 we can expect U E L"^ . 
Through a standard Fourier splitting argument we can only hope 

lk||2 <C(l+t)-3/4 

where the hold up for faster decay is the initial data. To get around this problem 
we will measure the difference 

w = u - $ where $ = e"^*/ 
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Here $ is the solution to the heat equation with initial data /. The function w 
satisfies a parabolic equation with zero initial data and a forcing term which we 
can control by restricting /: 

(2.4) wt + U -Vw + Vp ^ vAw - [/ • V$ 

W-w = Q w{0) = 

We can expect ||w(t)||2 to decay as (1 + i)^^/"* and if the heat flow corresponding 
to / decays at least as fast we can say the same about v, to make other parts of the 
argument work we need $ to decay faster. It is well known that the energy decay 
of the heat flow corresponding to / is intimately related to the behavior of / near 
the origin, therefore an assumption made on the decay of $ is really an assumption 
on / near the origin. With this in mind make the following assumption on /: 

Assumption 2.1. f (z X and there exists a po > such that /(^) = for every 
iei<Po. 

Remark 2.2. This assumption is really a bandpass filter for / which eliminates low 
frequencies. It is known that the corresponding heat energy decays exponentially, 
a fact demonstrated in the following lemma. Strictly speaking, one can relax the 
assumption on / so the heat energy decays at an algebraic (not exponential) rate 
and use the same method outlined in this paper, such / will need to be behave like 
a polynomial near the origin in Fourier space. See [1]. 

Lemma 2.3. /// satisfies Assumption \2.1\ and <& = e'^^^f, then 

(2.5) \ml<e-"'""'\\m 
Proof. The proof is quickly checked using the bound 

|$| = le-'-l^l'*/! 
<e--Pot|/| 

and computing the norm with the aid of the Plancherel theorem. □ 

3. Bounds for Stationary Solutions of the NSE 

Throughout this section we will assume / satisfies Assumption 12.11 and therefore 
$ = e^^^f satisfies ([13]). We are focused on the study of solutions for the two 
auxiliary PDEs: 

(3.1) W ■ + Vp = t/Af/'+i + / 

V • = 

and 

(3.2) wj+i + U' ■ Vw*+i + Vp = z/Aw^+i - W ■ V$ 

When dealing with either PDE we take the function W € fixed before hand. 
These PDE's will be used recursively to find approximate solutions for (|1.2p and 
p.4p respectively. In subsection 3.1 we recall existence theorems for these equations 
and Subsection 3.2 contains the decay rate calculations for In Subsection 3.3 

we make precise the notion [/* = v^it)dt which is then combined with decay 
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calculations in Subsection 3.4 to find uniform bounds on and show it is a Cauchy 
sequence in whose limit is a solution of (|1.2p . 

3.1. Existence Theorems. 

Theorem 3.1. Let W € and f G X . There exists a unique weak solution U^~^ 
to the PDE hS.l]) in the sense that for any € V, 

(3.3) < [/' • V[/'+\0 > +v < VC/'+\V(/) /,(/)> 
Moreover, this solution satisfies 

(3.4) l|VC/^+i||^< ^-^11/11^ 

Proof. We only outline the proof as similar PDEs are solved in the literature, see 
[5] , [8] , [11] , [16] , [19] , [20] , and [34] . A typical approach is to construct Galerkin ap- 
proximations {U^^}nGN by projecting the PDE onto finite dimensional subspaces 
of H^. A uniform bound similar to (|3.4p can be proved for each Galerkin approx- 
imation using an argument similar to that following (|2.ip . Once this bound is 
established it is possible to use the Banach-Alaoglu Theorem to find a subsequence 
{U^^}nGNcn that converges weakly in H^. The weak convergence is enough to 
pass to a limit in the linear terms. To pass through the nonlinear term one uses a 
stronger compactness theorem in the support of the test function 0. □ 

Theorem 3.2. Let U' G satisfy 

(3.5) ||VC/»|l2<^-^ll/llx 

and f satisfy Assumption \2. 1\ with $ = e'^^^f. There exists a unique weak solution 
^i+i g L°°(R+,L2) n L2(M+,7?1) to the PDE ( TOj) in the sense that for any 
0e Ci(R+;V), 

(3.6) < > + <U' ■ Vw'+\ > = -v < Vi(;^+\ V(/) > - < • V$, > 

V-w^+i^O w*+i(0)=0 
Moreover, this solution satisfies 

(3.7) sup||w'+i(t)||2 + W \\Vw^+\s)\\lds<Cpo^,^-'\\f\\j, 

t Jo 

Proof. The PDE in question is closely related to the Navier-Stokes equation and 
we refer to the literature for similar arguments, see [4], [5], [15], [16], [21], and [34] . 
It is typical to construct a sequence of Galerkin approximations which satisfies a 
uniform estimate similar to (j3.7p then use compactness arguments to pass through 
the limit. We give now a formal proof of p.7p which can be used as an a priori 
estimate in this approach. 

Multiply (|2.4p by w*"*"^ and integrate by parts, then use the bilinear relation (|2.2p 
to find 

- ||C/''||6||Vw;''+i||2||<f||3 
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The last line was obtained using Holder's inequality then Cauchy's inequality. 
Putting this together with the Gagliardo-Nirenberg-Sobolev inequality and the as- 
sumed bound on ||V?7'||2 yields 

(3.8) I||^^+i||2 + ,||v^^+i|l2<^||/||2^||<l>|l2 
Using Lemma [2?3l the bound (|2.5p implies 

Jo 2i^po 
Together with the Gagliardo-Nirenberg-Sobolev inequality and the heat property 
Io° \\^^i^)\\2 ds < ll/lll we estimate 

f\\m\\ids< f wmu^mhds 

Jo Jo 

< (^J\\m\\ldsy (^J\mis)\\lds 

(3.9) <Cp-K-'\\f\\\ 

Integrating (|3.8p in time then applying (|3.9p finishes the proof. □ 

Remark 3.3. In the theorems above the assumption W G is enough to ensure 
f/i . yu'+^ e (H^Y and W ■ Vw'+^ G (H^)' a.e. That is 

\<W ■WW+\<j)>\<C\\W(j)\\l and \ < W ■ Ww'+\(l) > \ < C\\W(j)\\l 

Therefore we are justified in multiplying the PDEs by t/*"*"^ and w'^'^^ respectively 
and integrating in space. Indeed, one just chooses a sequence of test functions ap- 
proximating either C/*+^ or and passes the limit through the weak formulation 
p.4p or (|3.6p . this will be a common technique in the remainder of the work. In 
both cases a stronger existence theorem is true but outside the scope of this paper. 



3.2. Decay of w. This subsection contains energy decay calculations for w, the 
estimates are an application of the Fourier Splitting Method with bootstrapping. 
The first step in the procedure is to apply the Fourier splitting method using the 
bound (|3.7p to find a preliminary decay rate. Once established, this preliminary rate 
is used to deduce a faster decay rate. This procedure is repeated until the recursion 
does not lower the rate again, in this case the hold up will be from estimates on 
the nonlinear term. The sequence of lemmas leading to Theorem 13.71 set up the 
bootstrap situation which is the main part of the proof for the theorem, establishing 
(I3.13P is the main goal of this subsection. The calculations are formal but can be 
made rigorous by applying them to a sequence of approximating solutions to (|2.4p 
(see [31]) or working directly with the weak formulation (see Remark 13. 3p . We 
begin with an estimate for \w\. 

Lemma 3.4. Let be the solution of US. ^|) given by Theorem \3.S\ with [/' and 

f satisfying the assumptions of the theorem. Then, 

(3.10) \w'+'\<C\mU'h(^l\\w^+\s)hds + ,.-'p^'\\f\\x^ 
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Proof. Through the Fourier transform of (13. 2p . noting the initial data is zero, we 
write 

i-t 

w'+^ = e-''l«l'(*-'')(e • U'w^+^ + + C • C/*$)(s) ds 
Jo 

Young's inequality with the Planchcrel theorem bounds 

+ \(M>\ < \\Wh{\\w'+^h + \m2) 
Taking the divergence of (j2.4p . then the Fourier transform, one can bound \p\ < 
C{\UHui+^\ + AU together, 



\w''''\<cmu''h / {w+'\\2 + \m2)is)ds 

Jo 

<CmU'h(^l\\w^+Hs)hds + p,',.-^f\\x^ 

The last line uses (|2.5p to evaluate the integral in time. This completes the proof. 

□ 

Lemma 3.5. Let w'"'^^ be the solution of US. g|) given by Theorem \3.2\ with C/' and 
f satisfying the assumptions of the theorem. Then, for any m > 4, satisfies 
the differential inequality 

l{ii+tr\\w^+Y2) 



(3.11) <cim,po,,^)\mm+tr-^i / \\w^+\s)hds+\\f\\ 



t X 2 

X 







c,.-^ffMi{i + ty 



Remark 3.6. In the statement of the lemma, the constant C(m,pojM) tends to c» 
as po — » or I' ^ and tends toward as ^ cxd. 

Proof. Multiply l\'d.2\i by it;*+^, after integration by parts then application of the 
bilinear relation (|2.2p and the assumed bound p.Sp we write 

(3.12) il||u;''+i||2 + j,||Vu;*+i||2 ^< . Vu;''+\ $ > 

< llC/HellVw^+i 



6||VW ||2||'?^||3 

iy3 



<-jll/llxll<i>ll3 + ;Tl|V^^ 



Now we split the viscous term in Fourier Space around the ball B{R) using the 
Plancherel theorem: 



JB(R)C 



B{B.)C 
JB(_R)C 

Jb{r) 
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Combining this with (|3.12|) : 



B(R) 



Then using p.lOp we bound 



B{R) 



<C\\W\\l[ f \W+\s)hds + u-'p^'\\f\\x] ( 



< C\ml { f \\w^+\s)hds + ^^-Vo ) R 







So, 



< C,yR'\\W\\l{l + v-'p^'f ( /* \W+\s)hds + WfWx 

+cu-^\\frx\mi 







In the preceding inequality we choose = —(1 + t)^^ then use (1 + i)™ as an 
integrating factor to estabhsh the lemma. Examining the line above one can see 
the constant in the statement of the lemma behaves like (1 + p^'^v~^)'^v~'^ /"^ , this 
is Remark 13.61 □ 



Theorem 3.7. Let he the solution of I2.4\l given by Theorem l 3. S\ with and 

f satisfying the assumptions of the theorem. Then, w*^^ satisfies the decay bound 

(3.13) \\uj'+\T)\\l < C(po,^)(l + + 11/111)11/111(1 + T)-^ 

Remark 3.8. The exponent of (1 + T) is such that ||w*'^^||2 is integrable over all 
time. The constant C{pQ, v) in the statement tends to oo as po ~* or ^ 0. It 
tends to as — + CX3 (see Remark [ 



Proof. Combining the bound on ||w'+^||2 given by (|3.7p with (|3.11|) we write 

I ((1 + trw+m) < dm, po, ,^)\miii + tr-'^ {t\f\\% + ii/iix) 

+ c(^)ii/iixii*ii3(i+<r 
< c(m,po,^^)(i + \m\i){i + 11/111)11/11^(1 +ir-^ 
+ c(:^)ii/iiiii'i>ii3(i+^r 
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The next step is to integrate in time, the first term on the RHS can be integrated 
directly while the second term is estimated similar to p.9p : 



<c{i^)\\f\\x r wmhwymui+trdt 




<C{po,,y)\\frx 
This gives an initial decay bound 

(3.14) \\w^+'{T)\\l < C{m,po, + ||C/1^)(1 + ||/||i)||/||^(l + T)-^ 

Now we begin the bootstrapping procedure. Proceeding in a nearly identical way 
to the argument immediately above, use (13. lip with ()3.14p instead of (13. 7|) . then 
integrate in time: 

||u;''+i(r)||^ < C(m,po,i^)(l + ||C/H^)'(1 + 11/11^)11/111(1 + T)-' 

This process can be repeated indefinitely but the "best" decay rate will be obtained 
after six iterations; here "best" is meant in the sense of best decay rate obtainable 
from p. lip . That this is in fact the best decay rate can be seen by examining the 
term ||u;*+i(s)||2 ds in (l3lT|) . once we have established ||w*+^||2 < C(l + i)~^ for 
fi > 1 this term integrates to a constant and we obtain the "best decay rate." As 
the bootstrapping steps are nearly identical to the above arguments and tedious to 
write out we skip to the final step: 

\\w^+\T)\\l < dm, po, ^.)(1 + imlfil + \\f\\x)\\.f\\xii + T)-^ 

□ 

Remark 3.9. In the above proof we relies on the exponential decay of $ which 
follows from Assumption l2.1l This can be relaxed to / ^ |^|'"' near the origin for a 
large k. 



3.3. Relation between (7* and w\ In this subsection we make precise, for our 
approximate solutions, the formal notion W = v'^(t) dt. We show approxima- 
tions of the integral v^{t) dt are bounded uniformly in and are Cauchy with a 
limit which is a solution of (|3.ip . Once this is established we apply the decay results 
from the previous subsection to find a uniform bound in for If. Throughout 
this subsection we use $ — e^*/- 

Lemma 3.10. Let be the solution of {2.4\l given by Theorem 1 3. '2\ with W and 

f satisfying the assumptions of the theorem. The function v^^^ — + $ satisfies 
!^v^+\t)dt€Ll. 

Proof. For each fixed i define the sequence {y„'+^}„gN C by 

Vt'= / dt 
Jo 
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Since ti'+^(t) G a.e. the sequence {V",j+^} is well defined. Relying on Minkowski's 
inequality for integrals with Assumption l2 . 1 1 (through (|2.5p ) and ()3.13p the following 
bound shows how the sequence {1^'*'^} is bounded uniformly (for n) in L^: 



ik:+'ii2< / \\v"^'m2dt 

Jo 



<C(po,^)(l + ||t/l^)=^(l + 11/111)^11/11 



X 



Similarly, 

rn+l 

(3.15) ||K|l-Kj+'ll2< / \\v'^%dt 



Observing ()3.13p and the decay of $ implied by Assumption l2 . II we know the integral 
/o°° Ik*"^^ll2 is finite so the RHS of (|3T^ tends to zero as n ^ CXI . Following well 
known arguments to prove a contraction lemma we can quickly deduce {K^^^jneN 
is Cauchy in and has a limit which we label u*+^(t) dt. □ 

Remark 3.11. The above lemma also implies v"^^^ (t) dt is finite a.e. in M.^. 



Lemma 3.12. Let be the solution of 112. 4\ ) given by Theorem 1 3. 2\ with and 
f satisfying the assumptions of the theorem. The function u*''"^ = + $ satisfies 
J^v'+^{t)dt = W+^ . 

Proof. To prove this lemma we show 1;'+^ (t) dt is a weak solution for p.ip then 
use the uniqueness implied by Theorem 13.11 to conclude the desired result. Let 
{V^+^}„gN be as in the previous proof. 

In (13. 6p choose cj) to be any member of V (so that it is constant in time). Use 
the relation v'"'^^ — + $ then integrate in time: 

^ < v'+\t),<l>> + <W ■ Vf'+i(t),0 > j dt 



(3.16) 



' / < Vw*+i(O,V0 > dt 
Jo 



After changing the order of integration and evaluating the first integral this becomes 

< v'+\n), (j)> + <W- S/V^+\(l) > = -iy< Vy„*+\ V0 > + </>> 

Observe the first term on the LHS tends to zero as n — s- cx). This follows from 
the decay bound p.l3p which implies w^'^^ tends to zero on compact sets and 
a similar well known property for the heat equation. The strong convergence of 
{^~'~^} enough to pass the limit through the remaining terms. Indeed, if 



< w ■ v(i/^+i - r+i),0 > I < c||vc/i2||i^^+' - r+i||2||v0|| 



3 



As n — > cx) this tends to zero for each test function (j> Q V, hence is a weak 

solution of (j3.ip . The uniqueness implied by Theorem 13.11 finishes the proof of the 
lemma. □ 
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Lemma 3.13. Let C/*+^ be the solution of U.S^) given by Theorem \3.1\ with W and 
f satisfying the assumptions of the theorem. Then the function C/*+^ satisfies 

(3.17) < Cipo, z.)(l + + \\f\\x)\\f\\x 

Remark 3.14. The constant C{po, v) in the above theorem tends to cxd as po ^ 
or ^ 0. It tends to as — + oo (see proof of Lemma [3?5|) . 

Proof. Define w*"*"^ = ui'"*"^ + Just as in the proof of Lemma 13.101 combine 
Minkowski's inequality for integral with (|2.5p and (|3.13p . but this time use the 
relation from Lemma 13.121 

\\U'+%< / \\v'+\t)hdt 
Jo 

< / \\w'+\t)\\2dt+ / mt)\\2dt 

Jo Jo 
<C{po,,^)il + \\Ur2fil + \\f\\x)^\\.f\\x 

□ 



3.4. Convergence of t/\ The goal of this subsection is to find the limit of the 
approximating sequence W and show this is a solution of the steady state Navier- 
Stokes equation. Later in the subsection we make two assumptions on /, they are 
smallness assumptions and allow a contraction argument to show is Cauchy. The 
assumptions will depend on how big we will allow the norm of U and throttle 
II V?7||2 so that a product of the and iJ^ norms of U is small. We will label this 
maximum value of the norm M (our choice) and keep it fixed throughout the 
remainder of this section. 

Lemma 3.15. Let f/*+^ be the solution of given by Theorem \3.1\ with C/' and 
f satisfying the assumptions of the theorem. There exists a constant C{po,h', M) 
so that ifWfWx < C{po,v,M) and \\U%2 < M then \\U'+^\\2 < M. 

Remark 3.16. For fixed po and i' the constant C(po, M) in the above lemma tends 
to as Af cxD or M 0. For fixed M the constant tends to as po or 
— )■ and tends to oo as ^ oo. 

Proof. By setting the RHS of (|3.17p equal to Af^ and considering Z — \\f\W ^ 
variable the proof is reduced to finding roots of the polynomial 

r^2 , y ^ T ^ 

C{p^,v){l + NPf 

Here, C{pQ,v) is exactly as in (|3.17p . Since i > This polynomial always has a 
strictly positive root, in this case the root is exactly the constant in the statement 
of the lemma. Indeed, 



(3.18) -^ + f^<Vl^ ^ 



2 - V^(a^(1 + M2)3 

Remark 13.161 follows by examining the RHS of (|3.18p . □ 



Theorem 3.17. Let M > and f satisfy Assumption \2. 1\ There exists a constant 
C{po, y, M) such that if \\f\\x < C(po, Af) the following hold: 
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(i) The PDE has a weak solution U G (in the sense of U.3]) ). 

(ii) This solution satisfies \\U\\2 < M in addition to i2. 1\) 

(iii) This solution is unique among all solutions which satisfy \2.1\) and have a 
finite norm. 

Remark 3.18. For fixed po and v the constant C{po, i^, M) in the above lemma tends 
to as M — > oo or M 0. For fixed M the constant tends to as po — > or 
v and tends to oo as — > oo. The behavior of the constant as M — > oo with 
the bound (|2.ip impUes ||J7||2|| VC/||2 < C(j^). In the time dependent Navier-Stokes 
system it is well known that when this product of norms is small for initial data 
the solution will remain smooth and well behaved, this result fits into that regime. 

Proof Choose U° S so that \\U°\\2 < M and \\VU°_\\l < i^-^WfWx- To construct 
such a function one could fix / then take a solution U for (jl.2p . At this point the 
solution is not known to be unique or have finite norm but by using a suitable 
cut-off function in Fourier space (C/° = x*U where x is equal to zero inside a ball 
containing the origin and one elsewhere) it is possible to limit the norm while 
not increasing the norm. 

Starting with U^, solve (|3.1|) recursively using Theorem 13.11 to find a sequence 
{U^}iZo which satisfies ||Vf7*||2 < Lemma [3.151 proviedes the uniform 

bound ||?7*||2 < M and so its limit, if it exists, must also satisfy this bound. We 
will now show this sequence is Cauchy in and a limit does indeed exist. The 
difference r*+i = - W solves 

After multiplying this by integrating by parts and using the bilinear relation 

(|2.2p one can deduce 

vW'^Y'+^Wl =< Y' ■ VY'+\U > 

< n6l|vy*+i||2||c/||3 
<^l|vri2||L/||2||vc/||2 + ^||vy^+i||^ 

The above sequence relies on Holder's inequality, Cauchy's inequality, and the 
Gagliardo-Nirenberg-Sobolev inequality. It implies 

\\VY'+^\\l < Ciy-^ WWhWVWhWVY'Wl 

(3.19) <Ciy-Hl\\f\\x\\WY'\\l 

Note the multiplication by is justified since all C/* (and hence all y) are 

bounded in H^. Using this bound recursively one finds 



\VY'+^\\l < {Ciy-Hl\\f\\xy+'\\VY^\\l 
< {C,y-Hl\\f\\xy+'2iy-^f\\ 



The last step relies the uniform bound on ||V?7'||2 < z^ ll/lli:- If WfWx < ^jj 
where C is the same as the line immediately above then Y^ tends to zero in 
which implies C/* is Cauchy, call its limit U. Through this construction we can also 
be sure ||J7||2 < M. Using standards arguments this strong convergence is enough 
to pass the limit through p.ip and show [/ is a solution of (|1.2p . For completeness 
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we will demonstrate how to pass through the nonlinear term: 

\<U ■ VU, (/)>-< C/^ • VU'+^,(j) >\<I + II 

1= I <C/-V(?7-C/*+i),0> I 
// = I < {U -W) ■WU'+\<j>> I 
To show / — > use Holder's inequality: 

/< ||t/||3||V(;7-C/*+l)||2||0||6 

Since the norm of U and the norm of (j) are bounded, the strong convergence 
W ^ U in shows the RHS tends to zero. The term // is handled in a nearly 
identical way. 

It remains to establish that U is the unique solution of (|1.2p among all solutions 
which satisfy (|2.1[) and have finite norm. Let U be any other solution which 
satisfies (|2.ip and has a finite norm. The difference Y = U — U solves 

(3.20) U -VY + Y -VU + Vp = iyAY 

U and U are bounded in and we are allowed to multiply this equation by 
Y. Then, proceeding in the same way as the lines leading to (|3.19p . 

\\\/Y\\l<C,.-Hl\\f\\x\\\7Y\\l 

The assumption on / made earlier in this proof is enough to guarantee \\f\\x < 
/CM and implies the solution is unique. □ 

Remark 3.19. This is exactly Theorem 11.11 Examining (|3.20p it seems Y G 
might be enough to obtain uniqueness as the first term on the LHS would formally 
integrate to zero after multiplying by Y . Unfortunately, following techniques used 
in this paper, we are not able to multiply (|3.20p by Y unless we know also Y G L^.. 
Indeed, using Holder's inequality and the Gagliardo-Nirenberg-Sobolev inequality 
one can see C/ • VF e (iJ^ n Ll)' C (H^)' , but not U ■ WY e {Hi)'. 

The uniqueness in this theorem, with the exact same proof, could instead be 
stated "This solution is unique among all solutions which satisfy (|2.ip and have a 
finite norm" and using some other technique it may be possible to expand this 
uniqueness theorem further. 



4. Stability of Solutions 

An important property of steady state solutions for physical problems is stability: 
"If a steady state solution is perturbed will it return to the same solution?" In the 
setting of the Navier-Stokes equation we investigate the stability of a solution U 
for (|1.2p by considering a perturbation wq and examining the long term behavior of 
the solution of the Navier-Stokes equation (|l.ip with initial data u(0) = JJ + wq- In 
particular, one would like to know what conditions on wq will guarantee solutions 
of (jl.ip approach U as time becomes large. An equivalent problem, found by 
subtracting ()1.2p from the (II. ip . is to determine when solutions of the following 
PDE tend to 0: 

(4.1) wt + u • Vu) + w • VC/ + Vp = i/Aw 

V • w = w{Q) = wq 
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A through examination of stabihty for the steady state Navier-Stokes equation is 
currently outside the reach of modern techniques (even in bounded domains) but 
using energy techniques we can prove strong stabihty results for perturbations of 
finite energy under certain restraints on /. Following the literature we introduce 
the following notion of stability, commonly called nonlinear stability, and find con- 
ditions on / which guarantee this type of stability. 

Definition 4.1. We say a solution U G of ()1.2p corresponding to f £ X and 
given by Theorem II. II is nonlinearly stable if it satisfies the following: if wq € is 
a perturbation and u is the solution of the Navier-Stokes equation (jl.ip with initial 
data Wq + U which satisfies, for any T > 0, 

u e L°°(0, T; L^) n L^{0, T; 

then 

(i) For every e > there is a (5 > such that 

||wo||2 < S implies sup \\u{t) — U\\2 < e 
teK+ 

(ii) u{t) tends to U as time becomes large, that is 

lim \\u{t) - U\\2 = 

t — >oc 

To start we can multiply (|4.ip by w, then integrate by parts and use the bilinear 
relation (|2.2p to find a formal energy inequality 

ljtM2 + '^W^Mll = - <w-\7U,w> 
< ||w||6||f/||3||Vu;||2 

<q|;7||j||vc/|i||ivHli 

The last two lines above were obtained with a combination of Holder's inequality, 
Cauchy's inequality, and the Gagliardo-Nirenberg-Sobolev inequality. If / is chosen 
so that U is small in norm (C|lC/|l2^^||V[/||2^^ < the above inequality 

becomes 

(4.2) ^lh«ll2 + ^l|Vii;||^<0 

This differential inequality implies ||w||2 is bounded uniformly. As ?7 € we can 
say the same about u = w + U. In other words, if one considers a solution U given 
by Theorem 11.11 then any finite energy perturbation will stay close to U, this is 
condition (i) in Definition 14.11 The rest of this section will be spent proving the 
stronger statement, that all finite energy perturbations return to U. 

4.1. Existence Tiieorems. Here we state the existence theorems and properties 
of the two PDEs examined in this section. Proofs of these theorems can be found 
in literature and are omitted. 

Theorem 4.2. Given T > 0, initial data uq G L^, and a forcing function f G 
the PDE (GUP has a weak solution 

(4.3) ueL^{{{),T]-Ll)C^L\[(),T]-Hl) 
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in the whole space M.^ which satisfies, for any e V, 

<ut,4>> +iy < Vm, V(/) > + < u • Vm, >=< f,4>> 

Proof See g], [5], [H], [H], [21], and [34]. □ 

Theorem 4.3. Lef u satisfy j| ) and ?7 e i/^. There is a constant C{v) such 
that if ||C/||i/i < C{v) the PDE {4.1^ has a unique weak solution 



w e L°^(R+,Ll)nL^(R+,Hl) 

satisfying, for any (j) G V, 

<wt,(j)> +v < Vw, V4> > + < u ■ Vw, (j)> + <w ■ VU, (p >=< /, ^ > 
as well as the following energy inequalities: 

(4.4) ^Ikll2 + '^I|V^«||^<0 

(4.5) snp\\w{t)\\l + ,, \\^w\\l<MO)\\l 
teK+ Jo 

Proof. The proof of this theorem is similar to the one immediately preceding with 
exception of the a priori bounds which are argued formally preceding (j4.2p . The lin- 
earity of the equation implies uniqueness (which is not known for general solutions 

of (in])). □ 



4.2. Decay of w. In this subsection we prove a decay property for solutions w given 
by Theorem 14.31 We show 11^112 ^ using a method developed for the Navier- 
Stokes equation in {25] . The method relies on generalized energy inequalities for 
the solutions which allow the energy to be decomposed in high and low frequencies. 
These are estimated independently and shown to approach zero. 

Lemma 4.4. Let u and U satisfy the assumptions of Theorem \4-.3\ with \\U\\h^ Isss 

then the given constant. Let cj) — e^'^' , ip — 1 — (j) and E{t) G C"'^([0, oo); L°°). 
The solution given by Theorem satisfies the following two generalized energy 
inequalities: 

U*w{t)\\l<\W''^^'-'^*w{s)\\l 

+ 2 I < M • Vw, e^''^^^-''^ * <^ * w > I dr 

(4.6) +2 I I < w • Vf7, e2''^(*~^)0 * * w > I dr 

J S 

E{tmw{t)\\i<muMs)\\i 

~2y [ E{T)U^w{T)\\ldT+ f E'{T)yMr)\\ldT 



(4.7) +2 J E{t)\ Kw-VU^ip^w > \dT 

+ 2 [ £:(t)| < iTW;, (1 - -02)i& > I dr 
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Proof. We give a formal proof here which can be made precise by considering an 
approximating sequence, see |25| for details. 

To see the first inequality multiply the PDE (|4.ip by e^''^'*"'"^^(^ * (j) * w and 
integrate from s to t. The assumptions are enough to ensure all integrals are finite 
and this multiplication makes sense. After integration by parts: 



U * < lle^^^*-").^ * w{s)\\l -V j II V(e''^(*-^),^) * w\\ldT 

+ 2 f \<u-Ww,e^''^^^-^'>^*^*w>\dT 

J S 

+ 2 [ I < w Vf/,e2''^(*~^)0*<^*w> Idr 



gi/A(t+s)^ describes a heat flow so the second and third terms on the RHS add to 
zero, this proves (|4.6p . For the second inequality, take the Fourier Transform of 
(|4.ip then multiply by tp'^w. After integration by parts one finds 



E{t)Uwmi<E{s)yjw{s)\\i-2iy / i?(r) || ^^^(t) || ^ dr 

J S 

+ I E'{T)\\%ljw{T)\\ldT + 2 I E{T)\<vj'^mj,'il:'^w>\dT 

J s J s 

+ 2 [ E{t)\ < -iTVw^tp^w > \dT 



The bilinear relation (|2.2p and the Planchcrcl theorem imply < u ■ Vw,w >= 
and (14. 7p follows immediately. □ 



Theorem 4.5. Let u and U satisfy the assumptions of Theorem \4.3\ with \\U\\hi 
less then the given constant and supj ||u||2 < oo. The energy of the solution given 
by Theorem\4.3\ decays to zero. That is, 



(4.8) lim||wm||2 = 

Proof. Following [25] we bound first the low frequencies using (|4.6p and then the 
high frequencies using (|4.7p and the Fourier Splitting Method. 

To show the low frequencies tend to zero we start by estimating the integrals on 
the RHS of 



I < It • Vw, e2''^(*~"')0 *<j)*w>\^\<<j)*<j)*u-w, e^''^^*-"') Vw > | 

< 110 * * U • U'||2|| Vu'||2 

< C||w||2||Vu;||2 

This estimate was obtained using integration by parts, the Cauchy-Schwartz in- 
equality. Young's inequality, and the Gagliardo-Nirenberg-Sobolev inequality. Sim- 
ilarly, 

\<w\/U, e'^''^^'-^^ * <^ * «; > I < CIIC/II2II V?«||^ 
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Combining these two bounds with (|4.6p yields: 



U * w{t)\\l < |le''^(*-^)<^ * w{s)\\l + C ( sup \\u{t)\\1 + \\U\\l ) f \\Ww\\l dr 

VreR+ / Js 

Heat energy is known to approach zero as time becomes large so 

limsup||(^*w(t)||^ <C( sup ||u(r)||^ + ||C/||^) [ ||Vw||^rfr 

The LHS is independent of s, noting the energy bound (|4.5p we see the RHS tends 
to zero as s — > cxd. Using the Plancherel theorem we conclude 

(4.9) lim \\(l3w{t)\\l = lim * w{t)\\l = 

t — ^oc t — >oo 

We begin work with the high frequencies on a similar path, bounding the integrals 
on the RHS of (gT]). Note ■0 = 1- e'l^l' £ L°°, then 

\ <w ■ VC/, i/'^w > \ = \ < S, - w - U, Tp'^w > I 

< W^hUwh 

< C\\w\\e\\Uh\\Vw\\2 
<C\\Uh\\^M\l 

This chain of inequalities used the Cauchy-Schwartz inequality, the Plancherel the- 
orem. Holder's inequality, and the Gagliardo-Nirenberg-Sobolev inequality. Simi- 
larly, but this time making use of the rapid decay properties of 1 — ■0^ , 

I < zTVw, (1 - 1p^)w > I < 11(1 - V^)w^||2||Cw||2 



< 



C||(l-^^ni6/5|k||2|lHl6||VHl. 



2 



<C||u||2||Vu;||2 
Use these two bounds with (|4.7p to find 
E{s 



(4.10) -l-cf sup \\u{T)\\l + \\U\\i) fE^WVwWldr 

Now split the viscous term and the term with E' around the ball with radius 
p(r) > 0, B{p): 



E{t)"^^ ^ J, E{t) 

1 

W) 



/ {E'{T)-2yE(T)p{Tf) I \{l~c^)w\'d(,dT 



B{p) 

^'^^^ f \{l-(j))w\^d^dT 
B{p) 



c 
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Upon choosing E{t) = (1 + and = a/2v{\ + t) (a> 3), (fiTTOj) becomes 

^11(1 

n/(\ 



1(1 -0Mt)||^< 1^^11(1 



(1 + 0" ys(p) 



+ C( sup h(T)||^ + ||[/||^) / ||Vu;|l^dT 

vreR+ 



Note |1 — 01 < I'CP if ICI < 1; so for large values of s, 



2 



|(l-0)z«|^de<C(l + T)-^||«;(T)||2 

S(p) 

This implies, again for large s, 

(4.11) < C sup ||u;(t)||^(1 + i)-2 

reR+ 

Taking into account the energy bound (|4.5p . this tends to zero as t becomes large. 
For any large s we are now justified in writing 

limsup||(l-(/.)u;(i)||2<cf sup ||u(t)||2 + ||[/||2') f \\S/w\\ldT 

Again relying on (|4.5p then letting s ^ oo we find 

hm \\{1 - (f>)wit)\\l = 

t — *oo 

Using this limit in the triangle inequality with the low frequency limit (|4.9p com- 
pletes the proof. □ 



Theorem 4.6. Let f satisfy the assumptions of Theorem \l.l\ and he such that ||/||x 
is less then the constant given by the theorem and \\U\\hi is less then the constant 
given by Theorem \4-.3\ The solution U of is nonlinearly stable in the sense of 

Definition\4.1\ 



Proof. Let u be given by Theorem l4.21 the difference v = u—U solves (|4.ip and u, 
U, f meet the criteria of Theorem l4.5l This proves (ii) in Definition l4.1l Integrating 
(|4.4p in time proves (i) in Definition 14.11 □ 

References 

[1] C. Bjorland and M. E. Schonbck. Poincare's inequality and diffusive evolution. Submitted, 
2007. 

[2] W. Borchers and T. Miyakawa. L^-decay for Navier-Stokes flows in unbounded domains, with 
application to exterior stationary flows. Arch. Rational Mech. Anal., If 8(3):273-295, f992. 

[3] W. Borchers and T. Miyakawa. On stability of exterior stationary Navier-Stokes flows. Acta 
Math., I74(2):3fI-382, 1995. 

[4] L. Caffarelli, R. Kohn, and L. Nirenberg. Partial regularity of suitable weak solutions of the 
Navier-Stokes equations. Comm. Pure Appl. Math., 35(6):771-831, 1982. 

[5] P. Constantin and C. Foias. Navier-Stokes equations. Chicago Lectures in Mathematics. Uni- 
versity of Chicago Press, Chicago, IL, 1988. 

[6] C. R. Doering and J. D. Gibbon. Applied analysis of the Navier-Stokes equations. Cambridge 
Texts in Applied Mathematics. Cambridge University Press, Cambridge, 1995. 



STEADY-STATE NAVIER-STOKES IN 



21 



P. G. Drazin and W. H. Reid. Hydrodynamic stability. Cambridge Mathematical Library. 
Cambridge University Press, Cambridge, second edition, 2004. With a foreword by John 

Miles. 

R. Finn. On steady-state solutions of the Navicr-Stokes partial differential equations. Arch. 
Rational Mech. Anal, 3:381-396 (1959), 1959. 

R. Finn. On the steady-state solutions of the Navier-Stokes equations III. Acta Math., 
105:197-244, 1961. 

R. Finn. On the exterior stationary problem for the Navier-Stokes equations, and associated 
perturbation problems. Arch. Rational Mech. Anal., 19:363—406, 1965. 

H. Fujita. On the existence and regularity of the steady-state solutions of the Navier-Stokes 
theorem. J. Fac. Sci. Univ. Tokyo Sect. I, 9:59-102 (1961), 1961. 

G. P. Galdi, J. G. Hcywood, and Y. Shibata. On the global existence and convergence to 
steady state of Navicr-Stokcs flow past an obstacle that is started from rest. Arch. Rational 
Mech. Anal., 138(4):307-318, 1997. 

G. P. Galdi and A. L. Silvestre. The steady motion of a Navier-Stokes liquid around a rigid 
body Arch. Ration. Mech. Anal, 184(3) :371-400, 2007. 

J. G. Heywood. On stationary solutions of the Navier-Stokes equations as limits of nonsta- 
tionary solutions. Arch. Rational Mech. Anal, 37:48-60, 1970. 

J. G. Heywood. The Navier-Stokes equations: on the existence, regularity and decay of 
solutions. Indiana Univ. Math. J., 29(5):639-681, 1980. 

E. Hopf. Ubcr die Anfangswcrtaufgabc fiir die liydrodyuainisclicn Grundglcichungen. Math. 
Nachr., 4:213-231, 1951. 

R. Kajikiya and T. Miyakawa. On decay of weak solutions of the Navier-Stokes equations 
in R". Math. Z., 192(1):135-148, 1986. 

T. Kato. Strong L^'-solutions of the Navier-Stokes equation in R"", with applications to weak 
solutions. Math. Z., 187(4):471-480, 1984. 

O. A. Ladyzhenskaya. The mathematical theory of viscous incompressible flow. Second Eng- 
lish edition, revised and enlarged. Translated from the Russian by Richard A. Silverman and 
John Chu. Mathematics and its Applications, Vol. 2. Gordon and Breach Science Publishers, 

New York, 1969. 

,1. Lcray. etude de divcrscs equations integrales non lincaires et de quelques problemes que 
pse riiydrodynamique. J. de Math. Pures et appl, 12:1-82, 1933. 

J. Leray. Sur le mouvement d'un liquide visqueux emplissant I'espace. Acta Math., 63(1):193- 
248, 1934. 

T. Miyakawa. On uniqueness of steady Navier-Stokes flows in an exterior domain. Adv. Math. 
Sci. Appl, 5(2):411-420, 1995. 

T. Miyakawa. On space-time decay properties of nonstationary incompressible Navier-Stokes 

flows in R". Funkcial. Ekvac, 43(3):541-557, 2000. 

T. Miyakawa and H. Sohr. On energy inequality, smoothness and large time behavior in for 
weak solutions of the Navier-Stokes equations in exterior domains. Math. Z., 199(4):455-478, 

1988. 

T. Ogawa, S. V. Rajopadhye, and M. E. Schonbek. Energy decay for a weak solution of the 
Navier-Stokes equation with slowly varying external forces. J. Funct. Anal., 144(2) :325-358, 
1997. 

C. W. Oseen. Neuere Methoden und Ergebnisse in der Hydrodynamik. Akademische Verlags- 
gesellschaft M.B.H., Leipzig, 1927. 

M. E. Schonbek. Large time behaviour of solutions to the Navier-Stokes equations. Comm. 

Partial Differential Equations, ll(7):733-763, 1986. 

M. E. Schonbek. Large time behaviour of solutions to the Navier-Stokes equations in 
spaces. Comm. Partial Differential Equations, 20(1-2):103— 117, 1995. 

M. E. Schonbek and T. P. Schonbek. On the boundedness and decay of moments of solutions 

to the Navier-Stokes equations. Adv. Differential Equations, 5(7-9) :861-898, 2000. 

M. E. Schonbek and M. Wiegner. On the decay of higher-order norms of the solutions of 

Navier-Stokes equations. Proc. Roy. Soc. Edinburgh Sect. A, 126(3) :677-685, 1996. 

M. E. Schonbek. decay for weak solutions of the Navier-Stokes equations. Arch. Rational 

Mech. Anal, 88(3):209-222, 1985. 

J. Serrin. On the stability of viscous fluid motions. Arch. Rational Mech. Anal, 3:1—13, 1959. 



22 



CLAYTON BJORLAND AND MARIA E. SCHONBEK 



[33] B. Straughan. The energy method, stability, and nonlinear convection, volume 91 of Applied 

Mathematical Sciences. Springer- Vcrlag, New York, 1992. 
[34] R. Temam. Navier-Stokes equations. AMS Chelsea Publishing, Providence, RI, 2001. Theory 

and numerical analysis. Reprint of the 1984 edition. 
[35] M. Wiegner. Decay results for weak solutions of the Navier-Stokes equations on R". J. London 

Math. Sac. (2), 35{2):303-313, 1987. 
[36] M. Wiegner. Decay estimates for strong solutions of the Navier-Stokes equations in exterior 

domains. Ann. Univ. Ferrara Sez. VII (N.S.), 46:61-79, 2000. Navier-Stokes equations and 

related nonlineaj: problems (Ferrara, 1999). 

Department of Mathematics, UC Santa Cruz, Santa Cruz, CA 95064, USA 

Department of Mathematics, UC Santa Cruz, Santa Cruz, CA 95064, USA 
E-mail address: cbjorlandQmath.ucsc.edu 
E-mail address: schonbek9math.ucsc.edu 



